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We show that strongly photoexcited graphene monolayers with 35fs pulses quasi-mstantaneously 
build up a broadband, inverted Dirac fermion population. Optical gain emerges and directly mani- 
fests itself via a negative conductivity at the near-infrared region for the first 200fs, where stimulated 
emission completely compensates absorption loss in the graphene layer. Our experiment-theory com- 
parison with two distinct electron and hole chemical potentials reproduce absorption saturation and 
gain at 40fs, revealing, particularly, the evolution of the transient state from a hot classical gas to 
a dense quantum fluid with increasing the photoexcitation. 

PACS numbers: 78.67. Wj, 73.22.Pr, 78.47. J-, 78.45. +h 



Graphene is gradually emerging as a prominent plat- 
form for ultrafast photonics and optoelectronics [l|-l3|. 
Growing evidence was demonstrated in, e.g., broad- 
band transparency and universal absorption from the 
near- infrared to visibleQ, carrier dynamics sat- 
urable absorption Q, pulsed photoluminescence 0-01, 
and coherently-driven photo-currents [lo| . For graphene 
to play a significant role in ultrafast laser technology 
or telecommunications that exceed semiconductor nanos- 
tructure performance, it is vital to investigate femtosec- 
ond nonlinearities of strongly photoexcited states. Prior 
time-resolved studies in graphene have been mostly con- 
cerned with the weak excitation regime where the pho- 
toexcited carrier density is much smaller than the initial 
background carrier density. In this case a linear power 
dependence of transient signals was observed 0, [lll - [l3| . 

Ultrafast photoexcitations strongly alter the thermo- 
dynamic equilibrium of electronic states and lead to 
a series of temporally overlapping rapid processes in 
graphene, as illustrated in Fig. 1(a). First, during or 
immediately following the pulse duration Tp, photoexci- 
tations are coherent. Then, electron-electron collisions 
lead to decoherence and eventually to a quasi-thermal 
transient distribution after a time rth- Finally, for longer 
times via coupling to phonons, the system relaxes back 
to equilibrium via cooling of the hot carriers (tc) and 
recombination of electron- hole pairs (tj.). In most semi- 
conductors and their nanostructures, where Tth S> Tp for 
~10 fs laser pulses jlj], ultrafast nonlinear photoexcita- 
tions lead to a largely non-thermal, peaked distribution 
close to the pump photon energy and state filling domi- 
nates on the 10s of fs time scale [first panel. Fig. 1(b)]. 

In this letter, we demonstrate that graphene is in the 
opposite limit Tth < Tp for 35 fs pulses, where a broad- 
band quasi-equilibrium, yet population inverted dense 
Fermi system emerges during the pulse propagation. The 
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FIG. 1. (color online) (a) Schematics of ultrafast optical 
interband excitations, (b) Dispersion of our electron-doped 
graphene monolayers {fi = 0.4eV) illustrating state filling 
(left) and band filling (right) that leads to stimulated emis- 
sion from a broadband, inverted population (red arrow). Also 
shown together is the pump pulse spectrum, (c) STM images 
of tomography of the sample used, (d) The static differential 
reflectivity spectra (red dots). The threshold AR/Rc for zero 
conductivity can be directly determined ~-1.46% (black line, 
see text for details). 



photoexcited phase space is quickly depleted via carrier- 
carrier collisions, leading to band filling [second panel 
of Fig. 1(b)]. Most intriguingly, we demonstrate that 
broadband gain emerges below the excitation photon en- 
ergy via stimulated emission. The broadband optical 
gain directly manifests itself via a remarkable negative 
conductivity at the near-infrared probe energy detuning 
even ^-^400 meV below the excitation level within hun- 
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dreds of fs. These results can be quantitatively described 
in terms of distinct electron and hole chemical potentials, 
which evolves from a hot classical gas to a dense quan- 
tum fluid with increasing the excitation from the linear to 
the highly nonlinear regime. Such femtosecond build-up 
of high-density and broadband population inversion has 
implications in advancing gaphene-based above-terahertz 
speed modulators, saturable absorbers and gain medi- 
ums. 

A Ti: Sapphire amplifier with central wavelength 
800nm {huj = 1.55eV) and Tp — 35fs is used to pump 
an optical parametric amplifier (900-2400 nm) to pro- 
duce near-infrared probes below the 1.55 eV excitation. 
Ultrafast degenerate and non-degenerate differential re- 
flectivity changes AR/R with ~ 40 fs time resolution are 
recorded with tunable pump fluences, which exhibits no 
pump polarization dependence. Our epitaxial graphene 
monolayer sample was grown by thermal annealing in ul- 
trahigh vacuum on a Si-terminated 6H-SiC(0001). The 
Fermi energy is ^ 0.4 eV in the sample (Fig. lb) reflect- 
ing the substrate-induced electron doping as reported 
15,[l3|. Fig. 1(c) shows the STM characterizations of to- 
mography of our samples used, which show homogenous 
carbon monolayer in atomic length scales (top) and 
scales with smooth overgrowth on the steps of SiC surface 
(bottom). The strong dependence of the 6a/3 x 6-\/3R30 
reconstruction modulation with the bias voltage con- 
firm the single layer graphene, with the homogeneity of 
the sample better than 90% across the entire probe re- 
gion. The monolayer thickness is further confirmed by 
the optical deferential reflection spectra [Fig. 1(d)], de- 
termined by the measurements with (R^+s) or without 
(Rs) graphene on SiC substrate. 

To understand optical response in graphene, we ex- 
pand the solution of the Fresnel equation with respect to 
a/c (of the order of the fine structure constant of quan- 
tum electrodynamics aqED = / {Tie) ~ 1/137 <C 1) 
and obtain (also see supplementary) 
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a' (w) =Re a{io) is the real part of the complex opti- 
cal conductivity. Eq.(IT]) leads to two key aspects: (i) 
the reflection coefficient provides a direct measurement 
a' (w) (or absorption) without reference to any model as- 
sumption,e.g., AR/R = 167r/(n^ - l)c* (cr' (r) - a' (0)); 
(ii) using the established universal value o'o = jf- for 
graphene monolayer without pulse [3, there exists 
a threshold value for the photoinduced differential re- 
flectivity that corresponds to the transition to a neg- 
ative optical conductivity and thus to gain behavior 
AR/R\^ = - ^"'^1° . With = 2.7 for SiC substrate 
follows AR/R\^ = —1.4582%., which can be determined 
experimentally by measuring differential reflectivity the 
optical spectra [Fig. 1(d)]. Here the reflection from 
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FIG. 2. (color online) (a) The degenerate differential reflectiv- 
ity at 1.55 eV for the graphene monolayer (solid dots and line) 
and SiC substract (empty dots) under the pump fluence of 
1842 /iJ/cm^. Inset: AR/R in logarithmic scale (black dots) 
together with the pump-probe autocorrelation (red line), (b), 
AR/R at 1.55 eV for different pump fluences. The dashed 
straight line is guide for the eyes. The different curves were 
taken in the order as marked. The perfect overlap of the pair 
of curves indicates the signals are not caused by laser-induced 
permanent changes, (c) The peak ^^\peak as function of the 
pump fluence (black squares) and the corresponding conduc- 
tivity change (red solid dots). Blue arrow marks the thresh- 
old for zero conductivity AR/R\^ = —1.4582% (see text). 
Dashed line: linear dependence (guide to the eyes). 



the zero conductivity state in the pumped garphene/SiC 
sample exactly corresponds to the case of bare SiC sub- 
strate. Critical to note also the negative a' (w) repre- 
sents the hallmark for the existence of gain in the ex- 
cited graphene layer, which should not be confused with 
the refection loss in the substrate. 

Next we present ultrafast degenerate reflectivity spec- 
troscopy to reveal femtosecond nonlinear saturation. A 
typical temporal proflle of AR/R at 1.55 eV in the 
graphene monolayer (black dots) is shown in Fig. 2(a), 
clearly showing a negative transient signal. Note the 
photoinduced change is negligible in the controlled ex- 
periment using the SiC substrate without graphene (red 
circles) . Several temporal regimes can be identified in the 
logarithmic scale plot in the inset: a pulse width limited 
rise ~ 40 fs (red line) , followed by a bi-exponential decay 
of 70 fs and 2.5 ps (blue lines). The power dependence 
of photoinduced AR/R at the maximum of AR in Fig. 
2(b) reveals a clear nonlinear behavior. Fig. 2(c) summa- 
rizes AR/R at the signal peak for different pump fluences 
Ip, showing a nonlinear pump fluence dependence above 
~ 1850^ J/cm^, at least one order of magnitude higher 
than what is reported for semiconducting single-walled 
carbon nanotubes 17[. Following Eq.([T]), the measured 



AR/R allows us to obtain the corresponding peak con- 
ductivity in photoexcited graphene, as shown by the red 
dot in Fig. 2(c) (normalized by (Tq)- At the highest pump 
fluence used ~ 3868/iJ/cm , AR/R peak appoaches 90% 
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of the critical value AR/R\^ and the peak conductivity 
drops to only 10% of (Tq- 

The most striking response is obtained after ultrafast 
non- degenerated differential reflectivity. Fig. 3(a) shows 
dynamics using 1.55 eV pump and low energy probes at 
1.16 eV and 1.33eV. It is clearly visible that, at high 
pump fluence, the critical value Ai?/i?|^, the threshold 
for negative conductivity (blue lines), indeed occurs for 
both cases. With increasing pump fluence. Fig. 3(b) 
indicates that, above Ipump,c — 2000/xJ/cm^, the stimu- 
lated infrared emission surpasses absorption loss in the 
photoexcited graphene for the 1.16 eV probe. The tem- 
poral profiles for pump fluence above the gain thresh- 
old indicate that the negative conductivity can persist 
for hundreds of fs, e.g., the 3960 ^J/cm^ at 1.16 eV. 
We emphasize three key aspects of this conclusion: (i) 
the critical value has not been reached by the degener- 
ate pump/probe [Fig. 2(c)] and appears exclusively for 
non-degenerate condition when probing below 1.55 eV, 
(ii) the AR/ R\^ is a model independent value that corre- 
sponds to (T = 0, which directly indicates the transition 
from loss to gain behavior, (iii) the femtosecond emer- 
gence of stimulated emission even at the ~400 meV below 
the excitation level indicates very rapid establishment of 
broadband population inversion and broadband gain in 
the strongly photoexcited graphene states. 

We can extract the number of photoexcited electrons 
in graphene immediately after the laser pulse from 
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where I{t,Ip) is the Gaussian pulse envelop I{t,Ip) = 
Ip\J exp ~4in2 ^2 ^ normahzed such that the total 
pulse fluence is Ip — ^I{t,Ip). For rth <C Tp, the 
absorption coefficient A (t) is determined by the adiabatic 
dependence of the absorption on the pump fluence and 
can be derived from the measured saturation curve at 
40 fs (Fig. 2c) by A(t) = + AA{t) = Ao(H- ^) 
(see supplementary). Without pump, the absorption of a 
graphene monolayer on SiC is A^ = (i+„3,c)2 Using 
the actual absorption A (t), instead of Aq, is crucial as the 
linear relation hajnex — Aolp substantially overestimates 
Hex during the pulse propagation, as shown in Fig. 3(c). 
From Eq.(3), we extract from the data extremely dense 
photo-excited fermions n^^P^ ~ 0.5 x 10^^ cm~^ for our 
electron doped sample, surpassing the saturation in a 
10 nm GaAs quantum well by more than two orders of 
magnitude under the similar excitation condition [l8|. 

Next we analyze the transient state at rth < t < Tc- 
Immediately after the pulse at At = 40fs, energy of the 
electronic system is conserved because no relaxation into 
the phonon system has taken place yet. In the case of 
highly excited graphene, the phase space constraint of the 
Dirac spectrum leads to an approximate conservation of 
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FIG. 3. (color online) (a) Ultrafast AR/R at 1.55 eV 
pump, 1.16 eV probe, at 1116 and 3960 /iJ/cm'^, and 1.33 
eV probe, at 4390 /iJ/cm^, respectively. Blue arrow marks 
AR/R\^ = —1.4582% for zero conductivity. Shown together 
are the pump and probe spectra, (b). The peak transient re- 
flectivity as function of the pump fluence. (c). The extracted 
transient fermion density at 40 fs (blue dots), as explained in 
the text, which is signiflcantly lower than Aolp/huj obtained 
from the universal conductivity (open circles), as illustrated 
in shadow area, (d). Theory (lines) vs. experimental values 
(rectangles) for non-degenerate (red) and degenerate (blue) 
transient conductivity at 40 fs. Shown together (lines) are two 
model simulations with the single (green) or distinct (black) 
chemical potentials. 



numbers of photoexcited holes and electrons, valid to the 
second order in the electron-electron Coulomb interaction 



lol . |20| . A recent explicit analysis of the short time dy- 
namics by Winzer et al. indicates that the conservation 
of hole and electron numbers is a good approximation for 
the high excitation regime 121|. Consequently, this gives 
rise to a slow imbalance relaxation and thus to a popula- 
tion inverted transient state with quasi-conserved occu- 
pations of the two branches of the Dirac cone in our ex- 
perimental condition. These, together with the assump- 
tion that a decohered, quasi-thermal state leads to the 
non-equilibrium transient distribution function {e = vp): 



fe(h) (e) - 
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characterized by the electron temperature, Te and two 
distinct chemical potentials /Xe and /x/i, for electrons in 
the upper and holes in the lower branch of the spectrum, 
respectively. Note that a scenario based on a single chem- 
ical potential does not explain the demonstrated popu- 
lation inversion. In thermodynamic equilibrium holds 
fie = —f'h = f- ^-iid = T. In general, Tg and fie(h) 
are functions of given photon energy, tku, total density. 
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FIG. 4. (color online) (a) The calculated transient electron 
(/ie), hole (/ih) chemical potentials and transient temperature 
(inset) as function of photoexcited carrier density, (b) Tran- 
sient electron and hole distribution function at 40 fs are plot- 
ted for different pump fluence. The two intersection planes 
represent the occupation probability of an electron-hole pair 
at the given excitation photon energy is equal to unity, i.e., 



and density of photoexcited carriers, riex (see supplemen- 
tary section). The non-equilibrium optical conductivity 
is calculated in Keldysh formalism [22| and follows as 



Ah 



(1 - /e {huj/2) ~ h {nLu/2)) . (4) 



Fig. 3(d) compares the peak transient conductivity 
with the calculated conductivity a (uj) of Eq|4] as func- 
tion of riox for two probe photon energies 1.55 eV and 
1.16 eV. Excellent agreement between experiment and 
theory is found which demonstrates the faithful repre- 
sentation of the transient state at 40 fs by the distri- 
bution function, Eq|31 For the degenerate scheme, our 
theory (black dashed line) yields cr — and thus per- 
fect transparency for n^'^^ — 0.48 x lO^'' cm^^. Once 
the systems is driven into this regime, a balance be- 
tween stimulated emission and absorption will lead to 
a transprancy. For non-degenerate scheme by probing 
at 1.16 eV, our theory (black solid line) predicts a crit- 
ical density n\l^ = 0.32 x 10" cm for the transition 
from loss to gain. All of these results agree quantita- 
tively with the experimental values = 0.5 x 10^^ 
cm~^ and nj^^^ = 0.34 x 10^* cm~^, respectively (black 
arrows) (23| . In addition, the experiment-theory com- 
parison of the conductivity is shown in Fig. 3d 
for the distinct- (black lines) and the equal-chemical- 
potential model (green lines) at the probe photon energy 
Tiui = 1.55eV and 1.16eV, which clearly identifies the va- 
lidity of the distinct- /X model calculation. 

The detailed information of the transient state at 40 fs 
as a function of ricx is shown in Fig. 4(a) for the transient 
chemical potentials and carrier temperature (inset). For 
weak excitation riox ^ w^/w^ (but larger than the num- 
ber of thermally excited carriers before the pulse) we find 
ksTe huj while < 0, corresponding to a hot, dilute 
gas of classical holes. Increasing n^x changes the sign of 
the hole chemical potential {^h > 0) and eventually de- 
creases Te, i.e. as a function of pump fluence we enforce 



a crossover from a hot and dilute Maxwell-Boltzmann 
gas to a degenerate population inverted quantum sys- 
tem, with Uf. -I- /ift, measuring the degree of population 
inversion [24]. Most notably, the transient conductivity, 
a (w), of Eq.dH) changes sign if + tJ-h crosses Tiuj, with 
broadband optical gain due to population inversion for 
the entire region w < /ie -I- /x/j (a < 0), as illustrated 
between two intersection planes in Fig. 4(c). The sep- 
aration between the planes is shown to approach to the 
pump energy 1.55eV at high excitation density, consis- 
tent with our experiment. 

We have showed the existence of pronounced fem- 
tosecond population inversion and near-infrared gain in 
strongly photoexcited graphene monolayers. These re- 
sults clearly reveal the transient electron and hole poten- 
tials are separated on the time scale of 100s of fs. Our 
experimental-theory comparison explains well of the ab- 
sorption saturation and gain from the dense fermions, 
and reveals a crossover from a hot Maxwell-Boltzmann 
gas to a degenerate dense Fermi-Dirac system. 
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